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ABSTRACT 


An experimental investigation was conducted to examine 
convection heat transfer in a duct of rectangular cross-sec¬ 
tion, having both a straight and curved flow passage. One 
wall was kept at a constant temperature, with the opposite 
wall being adiabatic. Heating of the air flowing in the 
channel was accomplished through Joulean heating of a wall 
composed of Temsheet. 

The experiments were conducted for steady state in 
Laminar and Transition flows. Development and progagation of 
the Taylor-Gortler vortices was shown to enhance the heat 
transfer rate in the curved section as compared to the 
straight section. 


4 








_ i inr m—nT ■ - .-dttir -m .. — —.— - 4' « » 

TABLE OF CONTENTS 

I. INTRODUCTION . 12 

A. TAYLOR-GORTLER VORTICES . 12 

B. HISTORY . 16 

II. INTENT OF THE STUDY . 2 2 

III. EXPERIMENTAL WORK. 24 

A. DESCRIPTION OF THE APPARATUS. 24 

B. EXPERIMENTAL PROCEDURES . 34 

IV. PRESENTATION OF DATA . 41 

A. ANALYSIS . 41 

B. RESULTS . 45 

V. DISCUSSION AND CONCLUSIONS . 51 

VI. RECOMMENDATIONS . 63 

APPENDIX A: EXPERIMENTAL UNCERTAINTY . 64 

APPENDIX B: SAMPLE CALCULATIONS . 69 

A. SAMPLE CALCULATIONS DATA . 69 

B. TEMPERATURE CALCULATIONS . 70 

C. POWER CALCULATIONS . 70 

D. AVERAGE HEAT TRANSFER 

COEFFICIENT (JT) . 72 

E. AVERAGE NUSSELT NUMBER (Su) . 72 

F. REYNOLDS NUMBER (Re) . 7 2 

APPENDIX C: CORRELATIONS . 73 

LIST OF REFERENCES . 76 

INITIAL DISTRIBUTION LIST . 80 


5 




























LIST OF TABLES 


TABLE I: SUMMARY OF STRAIGHT TEST SECTION RESULTS . 49 

TABLE II: SUMMARY OF CURVED TEST SECTION RESULTS . 50 

TABLE III: CORRELATIONS ..75 


6 











LIST OF FIGURES 


Figure 1 Schematic of Taylor vortices between 

cylinders ... 

Figure 2. Schematic of Taylor-Gortler vortices 

between parallel plates . 

Figure 3. Illustration of Equipment and Apparatus 

Figure 4. Cross-Sectional view of Channel . 

Figure 5. Straight Test Section, detailed schemat 

Figure 6. Curved Test Section, detailed schematic 

Figure 7. Thermocouple Placement in Test Sections 

Figure 8. Photograph of Test Channel . 

Figure 9. Straight vs. curved section results for 
the present study .. 

Figure 10. Comparison of present data with Durao, 
Ballard and McCuen for Laminar Flows, 
Straight Section .. 

Figure 11. Comparison of present data with Durao 
and Ballard for Laminar Flow, Curved 
Section.-.. 

Figure 12. Comparison of present data for Transitic 
Flows, Straight Section... 

Figure 13. Comparison of present data for Transitic 
Flows, Curved Section - 

Figure 14. Energy balance in straight section . 




















Symbol 


wo-wi 


K air 


Q air 


TABLE OF SYMBOLS 


Meaning 


cross sectional area of the channel 

area of the wall heater (Temsheet) 

specific heat of air at constant 
pressure 

height of the channel 

Dean number 

hydraulic diameter 

radiation shape factor 

average heat transfer coefficient 

thermal conductivity of air 

thermal conductivity of insulation 

mass flow rate of air 

local Nusselt number 

average Nusselt number 

average Nusselt number in curved section 

average Nusselt number in straight section 

Pradtl number 

atmospheric pressure 

rotameter flow percentage 

volumetric flow rate of air 

heat convected to the air 

heat lost through inner wall (Plexiglass) 

heat lost through outer wall (Temsheet) 

power supplied 


Units 


KJ/Kg °C 


2 °C 


W/m °C 
W/m °C 
Kg/sec 


3 / 

m /sec 






bymbol Meaning 

Q r heat transferred by radiation 

Q • heat transferred from the inner wall 

wi 

Q wo heat transferred from the outer wall 
R gas constant for air 

Re^ Reynolds number based on channel height 

Ri radius of curvature of inner wall 

Rpn electrical resistance of precision 

resistor 

R r total radiation resistance 

Ta Taylor number 


EXIT 


IN 


INS 


ROOM 


OUT 


T . 
wi 


wo 


PR 

H 


wi 


'wo 


bulk temperature of the flow 

flow exit temperature 

average flow inlet temperature 
temperature of the insulation 

ambient temperature 

average flow outlet temperature after 
each test section 

average temperature of inner wall 

average temperature of outer wall 

voltage across precision resistor 

voltage across wall heater (Temsneet) 

dimensionless axial length coordinate 

emissivity of inner wail (Plexiglas) 

emissivity of outer wall (Temsneet) 


Units 

W 

W 

w 

J/Kg ~K 


m 


R 


m 


-2 


°C 

°C 

°c 

°c 

°c 

u c 

°K 

°K 

V 

V 


9 


Symbol 

0 


Meaning 


Units 





f 

f 


! 

r 

h: 



b 


Y * 


air 


P 

a 


AT 

AT INS 

ax ins 


temperature correction to standard conditions 


dynamic viscosity of air Kg/m-sec 

densitv of air Kg/m 3 

2 4 

Stefan ? Boltzmann constant W/m °K 

mean temperature difference °C 

temperature difference in insulation °C 

thickness of insulation layers m 

2 

kinematic viscosity of air m /s 




i 

i 


{ 



10 







ACKNOWLEDGEMENT 


The author wishes to express his sincere appreciation and 
gratitude to Professor Matthew D. Kelleher for his invaluable 
assistance and direction in the formulation and development 
of this paper. Without his guidance this thesis would have 
been a near insurmountable effort. 

Additionally, the shop personnel deserve a special thanks 
for their assistance and encouragement, and without whom the 
physical set-up of the experimentation could not have been 
accomplished. 

This author must also pay a special thanks to his chil¬ 
dren, Keith and Robert, for their understanding in not having 
him home the many nights and weekends during the two years 
that this thesis culminates. And finally, to my wife, Susan, 
whose understanding, encouragement and love has been a 
source of inspiration and hope during this entire course of 
study. To her I must pay a very special thanks, she as truly 
has borne the burden of this academic endeavor and, without 
whom, this course of study and thesis could not have been 
accomplished. 


11 




I. INTRODUCTION 


A. TAYLOR-GORTLER VORTICES 

Since the early part of this century, considerable data 
and numerous studies have demonstrated that fully developed 
laminar flow along a concave wall does not remain two-dimen¬ 
sional [Refs. 1, 2, 3], Instead, the flow forms a system of 
spiral vortices, of counter rotating pairs, whose axes are 
aligned in the direction of the principle fluid flow. This 
phenomenon, known as Taylor-Gortler vortices, is the result 
of the variations in the centifugal forces acting on the 
fluid particles. Schematics of this type of motion are shown 
in Figures 1 and 2. 

In a channel that is curved in the streamwise direction, 
the fluid particles near the center of the flow cross-section 
are subjected to higher centifugal forces than those slower 
moving particles nearer to the boundary wall. Consequently, 
the tendency is for the fluid in or near the center of the 
channel to move outwardly towards the concave wall. This 
action then causes the fluid particles near the boundary wall 
to move in the spanwise direction, and finally move radially 
inward, replacing the outwardly moving particles. Once in the 
center of the channel these fluid particles come under those 
same stronger forces of the free stream velocity and continue 
the rotational effect. This cyclic motion causes the 
















Schematic of Taylor vortices between cylinders 
















formation, and propagation, of th.e counter-rotating Taylor- 
Gortler vortices-. 

Although the Taylor-Gortler vortices are primarily a 
laminar flow phenomenon they may affect the transition from 
laminar to turbulent flows, and are therefore important to 
understand {Ref. 2]. The additional mixing provided by the 
secondary motion has been thought to account for the increase 
in the rate of heat transfer from concave curved walls versus 
straight walls [Ref. 4]. 

The Taylor-Gortler vortices appear to have similarities 
to other vortex flow patterns such as the longitudinal vor¬ 
tex rolls developed in the laminar forced convection heating 
of fluid layers between parallel plates [Ref. 5]. It has 
also been thought that the cross-hatching, observed in 
reentry vehicles, can be explained at least in part, by the 
presence of streamwise vortices [Ref. &] . 

There are many possible applications that could result 
from a more thorough understanding of the Taylor-Gortler 
vortices and their effect on heat transfer and fluid flow 
characteristics. Applications, such as improved turbine 
blade cooling [Refs. 7, 8, 9], heat exchanger designs, as well 
as other such engineering designs and capabilities could take 
advantage of these vortices and improved heat transfer 
characteristics of the flow in the presence of these vortices. 








Man 


B. HISTORY 

The instability of an inviscid fluid rotating symmetri¬ 
cally about an axis, as in curved flow, was first considered 
by Lord Rayleigh in 1916 [Ref. 10], In assuming that the 
fluid was non-viscous, he determined that a necessary 
condition, for the stability of an inviscid fluid motion, was 
that the circulation increased with increasing radius. G. I. 
Taylor, in 1923, [Refs. 1, 11] expanded the earlier works of 
Rayleigh by an extensive analytical and experimental study of 
viscous fluids. His investigations focused on the flow 
between two cylinders, in which the inner cylinder rotated 
while the outer cylinder remained stationary. Taylor 
ascertained that such couette flows became unstable when the 
value of the dimensionless Taylor number, defined as: 



exceeded a critical value of 41.3. In the above equation, 
'd' is the width of the gap, assumed small when compared to 
'R^', the radius of the inner cylinder, ’U^' is the periphal 
velocity of the inner cylinder, and 'Re' is the Reynolds 
number. Taylor determined that for those cases in which the 
value of the Taylor number became greater than the critical 
value, a secondary motion developes and the Taylor vortices 
form. 
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Instability of a similiar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
considered analytically first by W. R. Dean [Ref. 12J in 
1928, for a channel formed by two concentric cylinders, 
where the radius of the inner cylinder was large in comparison 
to the small spacing between the inner and outer cylinder 
walls. Dean, in his studies, concluded that there would be 
an initiation of the flow instability, and the formation of 
vortices, (similiar to the Taylor vortices), when the Dean 
number, defined as: 

De = Re 

achieved a value greater than 36. In the above equation, the 
Reynolds number is based on the mean velocity of the unper¬ 
turbed flow, 'd' represents the channel half-width, and 'R^' 
is again the inner cylinder radius. The analytical work of 
Dean, was later verified by W. H. Reid [Ref. 13], using an 
approximate numerical solution. 

H. Gortler [Ref. 2], in 1940, studied the stability of 
laminar boundary layer profiles, on curved walls under the 
influence of small disturbances. In his study, Gortler found 
that these disturbances were similiar in nature to those 
studied by G. I. Taylor, which had led to the fluid flow 
instability in the form of vortices. By approximate numerical 
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calculations, Gortler concluded that only on the concave 
boundary walls were the amplified disturbances produced, and 
that the overall flow profile appeared to remain laminar in 
nature. In 1955, Gortler's approximate results were verified, 
with an exact solution, by G. Hammerlin, as reported by H. 
Schlichting [Ref. 14] and then further substantiated in a 
more extensive numerical analysis by A. M. 0. Smith [Ref. 3]. 
More recently, verification of these numerical solutions have 
been demonstrated by use of hot wire anemometry, laser doppler 
systems, and flow visualization techniques, [Refs. 15, 16], 
and in 1976, Y. Aihara [Ref. 17] conducted a non-linear 
analysis of the Gortler vortices. 

With the growing interest in the effects of the secondary 
flows associated with the Taylor-Gortler vortices, numerous 
specific studies have been published concerning the influences 
of these vortices on the transfer of heat in both the laminar 
and turbulent flow regimes. F. Kreith [Ref. 4], in 1955, 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent flows, 
and concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and under 
similiar turbulent flow conditions. 

In 1965, L. Persen [Ref. 18], considering the special 
cases of verv high and very low Prandtl number fluids, 
related the increase in heat transfer rate from a curved wall 





to the presence of the Taylor-Gortler vortices. His con¬ 
clusions were that the effect, of the development of the 
vortices, was to increase the rate of heat transfer through 
the boundary layer. 

One of the first experimental works on the effects of the 
Taylor-Gortler vortices with respect to enhanced heat transfer 
was by P. McCormack, et al., [Ref. 19], published in 1970. 

With the data, they concluded that it was necessary to retain 
the non-linear terms of the flow equation in order that the 
complete theoretical solution of the effects of the Taylor- 
Gortler vortices on heat transfer be explained. In 1977, 

R. Kahawita and R. Meroney [Ref. 20], furthered McCormack's 
work and concluded that the higher order terms, as well as 
the normal velocity components of the main flow, became 
increasingly more significant to the calculations at smaller 
wave numbers. 

In considering the forced convective heat transfer for 
fully turbulent flows, between horizonal parallel flat plates, 
Y. Mori and Y. Uchida [Ref. 5], in 1966, reported the 
formation of longitudinal vortex rolls with an axis parallel 
to the flow direction. These streamwise-axis directional 
rolls, similiar to the Taylor-Gortler vortices, formed when 
a critical value of temperature difference between the 
plates was exceeded. In 1971, M. Akiyama, et al., [Ref. 21], 
confirmed the formation of these vortices at a critical 
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Rayleigh number of 17Q8, for the case of heating from the 
bottom plate surface. 

There has been a rather limited amount of published 
literature regarding the flow and heat transfer in curved 
channnels of rectangular cross-section. Much of the pub¬ 
lished literature involves the development of numerical 
approximations and solutions for heat and mass transfer in 
curved ducts of various geometries (most notably, that of 
circular and triangular ducts). K. Cheng and M. Akiyama 
[Ref. 22] developed a numerical solution for forced 
convection heat transfer with laminar flows in curved channels 
of rectangular cross-section, but only foi small aspect 
ratios, (B/D; where, 'B' is the channel width and 'D' is the 
depth of fluid flow). Using the fully elliptical forms of 
the transport equations, G. Yee and J. A. C. Humphrey 
[Ref. 23], developed a numerical solution for flows and 
boundary conditions similiar to Cheng and Akiyama. K. Cheng, 
et al., [Ref. 24], in 1974, using numerical methods solved the 
Graetz problem for a curved channel of square cross-section 
CB/D-1). 

In 1963, B. S. Petukhov and V. N. Popov [Ref. 25], obtained 
an analytical expression for the Nusselt number for the flow 
of an incompressible fluid with variable properties in a 
circular tube. In 1976, A. A. Shibani and M. N. Ozisik 
[Ref. 26] using matched asymptotic expansion techniques for 
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the case of uniform wall temperature, solved the heat 
transfer problem between parallel plates with turbulent flow, 
for a wide range of Prandtl numbers. 

As to analytical and experimental studies of the parallel 
plate, channel flow heat transfer problem, Y. Mori, et al., 
[Ref. 27], obtained results for hydrodynamically fully 
de\eloped flows with constant wall heat flux, in curved 
channels of square cross-section. W. M. Kays and E. Y. 

Leung, [Ref. 28], experimentally obtained solutions for 
turbulent flow heat transfer in a concentric circular tube 
annulus with fully developed velocity profile and constant 
heat rate per unit length, for a fluid of Prandtl number 0.7. 
M. Durao [Ref. 29] and J. Ballard [Ref. 30], reported results 
for large aspect ratio (B/D = 40) rectangular cross-sectional 
channels with laminar flows. P. F. Brinich and R. W. Graham 
[Ref. 31], reported results for turbulent flows in a rectang¬ 
ular curved channel with an aspect ratio of 6, for the inner 
wall heated, the outer wall heated, and both walls heated. 

In 1979, F. D. Haynes and G. D. Ashton [Ref. 32], in studying 
heat transfer from a river to its ice cover, reported results 
for rectangular cross-section channels of aspect ratio 10, 
and compared them to the data obtained by Ashton [Ref. 33], 



IT. INTENT OF THE STUDY 


The primary purpose of this study was to investigate the 
effect of the Taylor-Gortler vortices, on the enhancement of 
th.e heat transfer rate in a curved channel of rectangular 
cross-section and large aspect ratio. This investigation 
began with flow velocities that were in the laminar flow 
regime and continued into the transition region, between 
laminar and turbulent flow. Further, this study compared the 
heat ransfer rate obtained in a curved section with the rate 
obtained in a straight section of identical aspect ratio. 

It was expected that the presence of the Taylor-Gortler 
vortices would enhance the heat transfer process at all 
velocities in the flow regimes investigated. The cause of 
this would be the secondary flow velocity components of the 
vortices, transporting heated fluid from the concave boundary 
wall inward, while simultaneously displacing the cooler fluid 
particles in the center, causing them to move outwardly 
towards the heated concave wall. 

This study was conducted using a single rectangular cross- 
sectional channel, that incorporated both a straight and a 
curved test section. The results obtained in the straight 
section were compared to those results obtained in the curved 
section, at the same flow velocities, in an effort to 
determine the effects of the Taylor-Gortler vortices on the 


f 

transfer of heat. Additionally, the results of both, the 
straight and curyed section tests were compared to the 
experimental results of Ballard JRef. 30] and Durao [Ref. 29]. 

Those results of the straight section that fell into the 
laminar flow region, were also compared to the analytical 
solution as given by P. McCuen, et al., [Ref. 35]; while 
those in the transition region were compared to the 
experimental results of Kays and Lueng [Ref. 28], and the 
analytical results of the Dittus-Boelter equation [Ref. 38] 
and Eckert's equation [Ref. 39]. While these equations are 
for fully turbulent flows in circular tubes they have been 
used to compare turbulent flows in other geometry channels 
when the diameter of the tube is replaced by the hydraulic 
diameter (wetted permimeter divided by the cross-sectional 
area). 

The straight section results, for transition flow, were 
also compared to the analytical solutions of Petukhov and 
Popov [Ref. 25] and Shibani and Ozisik [Ref. 26], the former 
for circular tubes while the latter is for parallel plates, 
both with fully turbulent flows. 

The curved section results for flows in the transition 
region were compared to the experimental results of Brinich 
and Graham [Ref. 31], and also to the solutions of Petukhov 
and Popov, and of Shibani and Ozisik, mentioned above. 





III. EXPERIMENTAL WORK 


A. DESCRIPTION OF THE APPARATUS 

A channel of rectangular cross-section built as described 
in references 29 and 30, was used to achieve the objectives 
of this study. The channel was manufactured from two 0.635 
centimeter thick sheets of plexiglas, separated by 0.635 
centimeter spacers whose inside lengths were used as the 
sides of the channel. 

The channel, as shown in Figure 3, was composed of a 
straight section, 122.0 centimeters long, followed by a 
curved section of 180 degree arc. The interior concave wall 
of the curved section has a radius of curvature of 30.5 
centimeters. The rectangular channel, with an aspect ratio 
of 40, was 0.635 centimeters high by 25.4 centimeters wide, 
with a cross-sectional area of 16.13 square centimeters and 
a wetted perimeter of 57.07 centimeters. A cross-sectional 
view of the channel, for both straight and curved sections, 
is shown in Figure 4. 

The working fluid used in this investigation was air, at 
room temperature, which entered the channel via an entrance 
bell constructed of Plexiglas, which was connected to the 
straight section. The entrance bell was designed and 
manufactured in accordance with ASME nozzle standards, with 
an elliptical curved base on a major axis equal to ten inches 
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and a minor axis of one inch. In an attempt to prevent the 
introduction of dust and other foreign matter into the flow 
channel, cheese-cloth was attached to the opening of the 
entrance bell nozzle. After passing through the channel, 
the flow was directed through an aluminum exhaust nozzle 
which was connected, by 3/4 inch PVC piping, to a Fisher and 
Porter Company variable area flow meter model 10A3565A. 

The rotameter, with a Bead Guide Tube number FP-1-35-G-10 
and a 1GSVGT float serial number 68T60, at 100% full scale 
had a flow rate of 0.793 cubic meters of air per minute 
(28.0 standard cubic feet of air per minute at one atmosphere 
and 70°F). 

The flow of air was drawn through the channel, piping and 
rotameter by a single electrically driven Spencer Turbo 
Compressor, 30 Hp. at 3500 rpm, rated at 550 cubic feet per 
minute at 70°F and one atmosphere. The flow rate was con¬ 
trolled by the use of valving both at the compressor and at 
the exit of the rotameter. 

To obtain the experimental heat transfer data for this 
study, two test sections of the channel were constructed: 

(1) a straight test section, 29.2 centimeters in length, 
located far enough downstream to ensure hydrodynamically 
fully developed flow, with a heated test section area of 
741.9 square centimeters. (2) a curved test section, 28.3 
centimeters in length, subtending an arc of 53.1 degrees, 
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located in the lower half of the concave curved portion of 
the channel, with a test section area of 717.4 square 
centimeters. 

In each of the test sections, the outer wall was 
modified as shown in Figures 5 and 6. Temsheet, a carbon 
impregnated porous paper with a property of uniform 
electrical resistivity, was glued to the interior surface of 
the outer wall of each of the two test sections. The flow of 
air was then to be heated by Joulean heating of the Temsheet. 
As the electrical resistance of Temsheet is not constant, but 
rather slowly variable with temperature, a precision resistor 
with an electrical resistance of 2.078 ohms was connected in 
series with the Temsheet to allow calculation of the 


instanteous power being supplied. 

The variables to be measured, in this investigation, were 
(1) the inlet flow temperature for each of the test sections 
(Tj^), (2) the outlet flow temperature at the exit of each of 
the test sections (Tq^^) > (.3) the temperature of the heated 
boundary wall of each of the test sections (T wo ) » (4) the 
temperature of the unheated boundary wall of each of the test 


sections (.T^^) , (.5) the flow temperature at the exit of the 
channel (T EXIT ), (6) the temperature between the two insula¬ 
tion layers for each of the two test sections (Tj^g), (7) 
the voltage across the precision resistor (V p ^), (8) the 
heater voltage (V^), and (9) the flow rate through the 
rotameter (Q). 
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^Temsheet 





6. Curved Test Section, detailed schematic 







The temperature measurements required ahoye, as well as 
all other temperature measurements in this investigation, 
were made with Copper-Constantan thermocouples. All the 
thermocouples had been previously calibrated using a 
ROSEMOUNT Communtating Bridge model 920A, and a ROSEMOUNT 
Constant Temperature Bath model 913A, as outlined by Ballard 
[Ref. 30]. 

A total of sixty-four Copper-Constantan, glass insulated, 
30 gauge thermocouples were constructed and inserted in 
specific locations to facilitate the measurement of the 
required temperatures. All the thermocouples were wired to a 
Hewlett Packard Data Acquisition System, model 2010c, which 
enabled an automatic printed record of the instantaneous 
reading for each thermocouple location. All the thermocou¬ 
ples were also connected to a Kaye Instrument Ice Point 
Reference, model K140-4, serial 1780. 

In each of the test sections, five of the thermocouples 
were inserted between the first and second layers of 
ARMAFLEX 22 insulation, and another five inserted between the 
second and third layers. Each of these sets of five 
thermocouples were connected in parallel to facilitate the 
recording of an average value of the temperature (TINS) 
between the respective layers of insulation. Additionally, 
sets of four thermocouples, also connected in parallel to 
indicate average temperature, were inserted into the channel 
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at the following three locations; l.i the entn i: ■' of the 
channel, (2) the exit of the siraigh: . ■■■- r section, and (5) 
the exit of the curved test section. At ; ach or' these three 
locations, the thermocouples were se.' five -entimeters 
apart in the spanwise direction to a Ll--w reading of the 
bulk temperature of the flowing air. 

In each of the two test section-, eight thermocouples 
were placed in direct contact with the Temsheet through 
small diameter holes drilled into the Fiexiglas outer walls, 
with the thermocouple beads electrically insulated from the 
Temsheet with ENMAR Heat Resisting Glyceryl Phthalate. The 
final sixteen thermocouples were placed, eight to a test 
section in the Plexiglas unheated boundary wall, opposite 
those positioned in the Temsheet. A detailed sketch of the 
thermocouple placement is provided in Figure 7. The average 
heated outer wall temperature (T ), was provided by those 
thermocouples in contact with the Temsheet. The average 
unheated inner wall temperature (T w ^j, was provided oy tne 
tnermocoupies embedded in the Plexiglas wall. 

Tne thermal insulation consisisted of layers of 
ARMSTRONG ARMAFLEX 22 Sheet Insulation, a flexible foamed 
plastic material, on the heated side of each test section, 
and black foam insulation one-half inches thick elsewhere. 
Three layers of the ARMAFLEX 22 insulation were used, each 
layer had a thickness of 0.636 centimeters, and covered an 
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area slightly larger than the respective test section. 

These layers of insulation were used in the compution of 
heat flux losses to the environment, sample calculations 
of these losses are given in Appendix C. The rest of the 
channel was covered with the one-half inch thick Black 
foam insulation to reduce the losses to the surroundings 
through the plexiglas walls as the heated air transited 
the channel. The insulation was held in place by the use 
of Adhesive Heat Resistant ventilation Duct Tape as seen 
in Figure 8. The entire length of the channel as well as 
all the connections between the channel and the rotameter, 
were sealed with General Electric Silicone Rubber Sealant 
Caulk to ensure that there was no leakage of air into the 
channel or rotameter downstream of the entrance bell 
nozzle. 

Two aluminum electrodes, 0.318 centimeters thick, were 
inserted in between the Temsheet and the plexiglas, in each 
of the test sections, as shown in Figure 4. The plates, 
with the Temsheet and aluminum electrodes were assembled 
in the channel with plastic teflon screws for the purpose 
of electrical insulation. 

B. EXPERIMENTAL PROCEDURES 

Experiments were conducted with the volumeteric flow 
rate varied incrementally from 0.198 to 0.674 cubic meters 







Test Channel 









of air per minute to correspond to particular Reynolds 
numbers. The data was recorded and calculations performed 
for each flow rate. The experimental procedures followed 
were the same for both the straight ana curved test 
sections. The straight test section runs served as a 
baseline for the curved test section and allowed tor 
comparison of the results. The appropriate electrical 
power required to heat the Temsheet boundary wall to 
approximately 50°L was determined through preliminary 
experiments. This temperature was choosen to allow at 
least a twenty degree celisus temperature difference 
between the unheated and heated boundary wall average 
temperatures and to facilitate achieving steady state 
conditions within a reasonable time span. Estimates of 
time to achieve steady state conditions and decisions as 
to data collection were also based on these preliminary 
test runs. 

The determination as to what constituted a steady 
state condition required careful consideration in view 
of the multiple input variables and their effects on the 
results. The final criteria for steady state was based on 
three input variables: (1) flow rate percentage, as 
recorded from the rotameter; (2) heated boundary wall 
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temperature, as recorded from the thermocouples in that 
area; and (3) the heater voltage, as recorded on the 
voltmeter. When these conditions varied by less that 
two percent over a ten minute period, the experiment 
was considered to have achieved steady state. Under 
these conditions the assumption of steady flow and 
constant heat flux from the heated wall were assumed 
valid. 

Originally, the time required to reach this steady 
state condition was found to be in the range of two 
hours, accordingly all runs were at two and one-half 
hours durations. Data was taken on one hour intervals 
for the first two hours and then at ten minute intervals 
for the next thirty minutes. Later it was found that the 
higher flow percentages for the curved test section, 
required about thirty minutes more time to full achieve 
steady state, so test runs with a flow percentage of 
sixty percent or greater were extended with three 
additional ten minute interval readings for a total time 
of three hours. 

The instantaneous values for all the temperatures 
were automatically printed by the digital recorder in 
millivolts. A standard thermocouple table was used in 
conjunction with previously obtained thermocouple 
calibration data to convert the millivolts to degrees 
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Celsius. The volumetric flow rate was read directly from 
the installed rotameter as the percent of full scale flow. 
The instantaneous power supplied was calculated using the 
relationship: 

V H V PR 

Q p = - 

r pr 

Where the precision resistor voltage (V pR ) and the heater 
voltage (V^) were recorded in volts from the installed 
voltmeter in the digital Data Acquistion System, and the 
electrical resistance of the precision resistor (Rp p ) had 
been previously determined to be 2.078 ohms. 

As the area of the curved test section was slightly 
less than the area of the straight test section, (0.0718 
vs. 0.0742 square meters), the heat flux from the heated 
wall, of the respective test sections was different. The 
room temperature entering the test section was checked 
with the use of a mercury thermometer located in the 
vicinity of the apparatus. 

During the preliminary test runs of the straight 
section, the thermocouple array recording the exiting 
temperature of the flow of air, was set at three 
millimeters from the heated wall. This was satisfactory 
at the lower volumetric flow rates, but at the higher 
flow rates it led to calculations in which the power 
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convected to the air was greater, by approximately fifteen 
to twenty percent, than the power being supplied to the 
heater. It was determined that the thermocouples were 
not indicating the correct bulk temperature of the 
exiting flow of air, but rather a higher value of temp¬ 
erature, this then resulted in an erroneously high value 
of convected heat transfer, and higher Nusselt numbers. 

To correct this, the system was allowed to run until the 
walls of the downstream curved section achieved steady 
state. Then, calculating the losses, by use of the 
thermocouples in between the insulation layers of the 
curved test section, and applying these losses to the 
entire curved section area, the energy loss was calculated. 
The flow exiting the channel was assumed to be thoroughly 
mixed and therefore the temperature as recorded by the 
positioned thermocouples, at the exit of the curved 
section, was used as an accurate representation of the 
average bulk exit temperature. With the final exit 
temperature and the losses now known, the temperature of 
the air as it exited the straight section was calculated. 

To correct this condition, the thermocouples were 
inserted further in to the channel, or raised closer to 
the outer wall until the millivolt reading was compatible 
with the calculated correct temperature. After this 
procedure, accomplished at a flow rate of 0.595 cubic 
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centimeters of air per minute (751 of full scale) was 
completed, additional test runs were conducted at 0.555 
and 0.635 cubic centimeters of air per minute (70% and 
80% respectively), without adjusting the thermocouple 
positions, to check the accuracy. In each case, the 
temperature was within 0.2S o C of what the energy balance 
indicated it should be. 

It was later discovered that the curved section re¬ 
sults had the same problem at flow rates greater then 
0.595 cubic centimeters of air per minute. Since the 
channel construction is such that it does not allow 
further downstream data collection, and since the power 
difference was less than four percent it was decided to 
proceed with the correlation of the data and experimenta 
tion and note the discrepancies. 
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IV. PRESENTATION OF DATA 


A. ANALYSIS 

The uniform electrical resistivity of the Temsheet, 
approximated a constant heat flux surface for the heated 
wall. With insulation on the unheated wall the heat losses 
through that wall were negligible and the wall could be 
considered adiabatic. With the large aspect ratio of 40, 
the channel was considered to be approximately parallel 
plates. The straight portion of the channel upstream of 
the straight test section, was of sufficient length to 
ensure that the flow was hydrodynamically developed for 
the flow velocities investigated. Based on these assump¬ 
tions, the experimental configuration was modeled as 
forced convection between parallel plates, for hydro- 
dynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant specified heat flux, and the 
other wall insulated. 

To analyze this condition, several quanties were 
defined as follows: 

The average heat transfer coefficient between the 
heated wall and the flow of air was defined by the equation: 

Q . hA pL AT 
air 








f 


where, '0 . ' was the heat convected to the air, 'H - ' 
was the average heat transfer coefficient, 'Ap^' was the 
area of the heated Temsheet for the section being tested, 
and 'AT' was the difference between the average heated 
wall temperature (T ) and the average bulk temperature 
of the fluid (Tg). The average bulk temperature of the 
fluid was defined as the arithmetic mean of the fluid 
inlet bulk temperature (T-^) and the fluid outlet bulk 
temperature (Tq UT ). The delta temperature (aT) could then 
be expressed as: 

AT = T- T = T - Hi_ 

wo B wo 2 

The actual heat convected to the air (Q . > was 

air ) 

calculated separately using th? expression: 

Q - m c (T - T ) 

air ^ air OUT IN 

where ' c p a p r ' was the specific heat of the air at constant 
pressure and 'm' was the mass flow rate of the air. The 
mass flow rate, is further defined as the product of the 
volumetric flow rate of the air and the local flow density 
of the air as it transits the rotameter. The volumetric 
flow rate of the air was corrected to standard conditions 


for the rotameter and the local flow density was calculated 






assuming "perfect gas" behavior at constant pressure and 
utilizing the temperature of the air as it exited the 
channel CT^xiT)* 

The average Nusselt number was then calculated from the 
average heat transfer coefficient, with the expression: 

R d 

mr = - 

K air 

in which ’d’ is the height of the channel and 'K . ’is the 

air 

thermal conductivity of the air evaluated at the average 
bulk temperature. 

The Reynolds number was then calculated for each test 
run as follows: 


m d 



where again 'm' and ’d' are the mass flow rate and height 
of the channel respectively, 'A ' is the cross-sectional 
area of the channel, and ’ u' is the dynamic viscosity 

all 

of the air evaluated at the exit flow temperature. 

For the experiment runs utilizing the curved section, 


the Dean number, defined as: 




was also evaluated. In this equation, 'd' is again the 
height of the channel, and *R^ * is the radius of curvature 
of the unheated wall surface of the channel (the convex 
wall) . 

The heat losses, by conduction, through the heated 
plate and the insulation, by conduction, were also cal¬ 
culated for each test section utilizing the expression: 

AT 

INS 

Q = - 

lo ( x /K ) A pt 
INS INS 


With 'ATj^g' being the difference in temperature, as record¬ 
ed, oetween the first and second layers of AkMAILEa 22 
insulation, and that recorded between the second and third 
layers. f Xj^ s ’ was the thickness of each layer of insula¬ 
tion, * ^INS * was t * ie 'thermal conductivity of the 
insulation, and T Ap L ’ was again the area of the heated 
surface. 


In order to verify that Q„. , as calculated bv the 
temperature rise in the air, was actually the convection 
heat transfer from the Temsheet heated surface, an 
estimate of the radiation heat transfer was also calculated 
for each test run. The eight thermocouples in the lower 
Plexiglas surface measured the temperature of that surface 
and allowed the computation, utilizing the Stefan-Boltzman 
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law of radiation, of the radiated heat transfer (Q ) with 
the expression: 


^r = 


ff (T 


wo 


T wi> 


R 


R 


where 'cr * is the Stefen-Boltzman constant, and 'R * is the 
total resistance to radiation heat transfer between the 
surfaces as computed by: 



^PL e wo fe wi 


In the above equation, ’« wo ’ is the emissivity of the 

heated Temsheet wall surface and 'e . ' is the emissivity of 

wi 

the unheated adiabatic wall surface (Plexiglas), 

A sketch of the control volume and the energy balance, 
with sample calculations for one test run, is given in 
Appendix B. 


B. RESULTS 

The data obtained from the experimental test runs were 
evaluated utilizing the expressions described in the 
previous section, and as illustrated in Appendix B. The 
results of the major variables are shown in tabular form, 
in Table I for the straight test section, and Table II for 
the curved test section. A plot of the average Nusselt 
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number verses Reynolds number, is given in Figure 9 for 
comparison of the two test sections. The uncertainty bands 
have been indicated for the plotted experimental data and 
a sample calculation for the complete error analysis is 
given in Appendix A. 

As shown, the results indicate an increase in the 
rate of heat transfer with increasing Reynolds number for 
both the straight and curved test sections. Additionally, 
the heat transfer rate of the curved test section was 
higher than that of the straight test section for each 
Reynolds number investigated. 

In investigating the results obtained, and considering 
those Reynolds numbers below 2000 (1025 based on this 
channel configuration) as laminar, at least squares 
straight line correlation of the data results in: 

0.67 

Flu = 0.06 5 Re d 

for laminar flows, and 

— 0.58 

Nu = 0.117 Re d 

for transition regime flows, in the curved test section. 
With correlations of 

__ 0.25 

Nu = 0.983 Re ^ 
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for laminar flows, and 


— 0.56 

Nu = 0.110 Re d 

for the transition flows in the straight test section. 

In the laminar flow regime, the increase in heat 

transfer was approximately fifteen percent, while in the 

higher transition flow region the increase was approximately 

thirty percent. This compares favorably to the earlier 

studies by Kreith [Ref. 4], in which he reported that the 

heat transfer rate along concave walls increased from 

twenty-five percent to sixty percent for Reynolds numbers, 

4 6 

based on hydraulic diameter, between 10 and 10 , and 
Ballard [Ref. 30], who found an eleven percent increase 
in heat transfer rate for laminar flows. 
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Figure 9. Straight vs. Curved Section results for the Present Study 






TABLE I 




SUMMARY OF STRAIGHT TEST SECTION RESULTS 


r 

\ 

Re d 

Qair 

(W) 

IT 

(W/m 2 °C) 

AT 

(°C) 

>Hi 

i 

[f. 

784 

39.42 

21.28 

26.7 

5.08 

' 

941 

40.14 

22.56 

24.0 

5.38 

i 

1099 

43.26 

23.11 

25.2 

5.51 


1257 

45.13 

23.93 

25.4 

5.71 

•4 

1414 

51.80 

25.89 

27.0 

6.18 

‘1 

1571 

50.44 

26.59 

25.6 

6.34 

- 

1729 

53.50 

28.61 

25.2 

6.83 


1886 

57.01 

31.32 

24.5 

7.47 

► 

2042 

59.24 

31.99 

25.0 

7.63 

« 

t 

2203 

61.65 

33.50 

24.8 

7.99 

■ 

2361 

62.97 

35.30 

24.0 

8.42 


2517 

67.11 

35.44 

25.5 

8.45 


2675 

69.14 

36.39 

25.6 

8.68 
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TABLE II 




SUMMARY OF CURVED TEST SECTION RESULTS 


Re d 

De. 

Qair 

(W) 

K 

(W/m 2 °C) 

AT 

(°C) 

Nu 

782 

114 

41.19 

23.63 

24.3 

5.64 

937 

137 

47.18 

26.57 

24.7 

6.34 

1092 

159 

50.88 

28.56 

24.8 

6.81 

1172 

171 

53.71 

30.49 

24.5 

7.27 

1251 

183 

57.16 

32.32 

24.6 

7.71 


i 1406 

205 

58.11 

33.01 

24.5 

7.88 

1570 

229 

62.63 

35.12 

24.9 

8.38 

1720 

251 

69.41 

37.23 

26.0 

8.88 

1883 

275 

70.90 

38.29 

25.8 

9.13 

2038 

297 

76.55 

41.66 

25.6 

9.94 

2196 

320 

80.36 

42.72 

26.2 

10.19 

2351 

343 

83.30 

45.98 

25.2 

10.97 

1 2505 

366 

88.82 

46.87 

26.4 

11,18 

2659 

388 

94.67 

49.88 

26.4 

11.90 





V. DISCUSSION AND CONCLUSIONS 

With the measurement of the unheated wall and the air flow 
temperatures, it was indicated that there was a negligible 
difference between the fluid bulk temperature and the tempera¬ 
ture of the unheated wall. Additionally, the energy balance 
and sample calculations, provided in Appendix B, demonstrate 
that the radiated heat transfer from the unheated wall sur¬ 
face is minimal. Thus it was assumed that the heat transfer 
to the flowing air was solely by convection from the heated 
Temsheet. 

The high aspect ratio of the channel as discussed pre¬ 
viously, provided the original basis for the assumption that 
the experimental apparatus, as configured, was one of infinite 
parallel plates. The experimental data, for those Reynolds 
numbers in the laminar flow region, tend to substantiate this 
assumption, in that they approach the theoretical limit for 
average Nusselt number of 5.385, for parallel plates with one 
wall heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 37]. For a Reynolds number, based on hydraulic 
diameter, of 1835, the average value of Nusselt number obtained 
was 5.382. 

With the formulation of the problem complete, the 
experimental data obtained can now be compared to the analytical 
solutions and experimental results of the same or similar 
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problem. That problem being, flow between infinite parallel 
plates with one wall at constant heat flux and the opposite 
wall adiabatic. 

A comparison of the experimental results obtained in this 
study and the experimental results of Durao [Ref. 29] and 
Ballard [Ref. 30] (who investigated the problem solely for 
laminar flows), are shown in Figure 10, for the straight sec¬ 
tion, and Figure 11, for the curved section. In both these 
studies, as in this investigation, the channel has an aspect 
ratio of 40, and the experimental procedures were similar. 

It can be seen that the results obtained in this study, in 
the laminar flow regime, appear to be in agreement with the 
data obtained previously for both test sections. Additionally, 
numerical correlations of the present data, as given earlier, 
for the laminar flows compares favorably to the numerical 
correlations of both Durao and Ballard, for Reynolds numbers 
greater than 600, as given in Appendix C. In the two investi¬ 
gations, mentioned above, comparisons were made to the analy¬ 
tical studies of McCuen, et. al. [Ref. 35], for heat transfer 
between infinite parallel planes with constant wall temperatures 
and heat flux. Ballard also compared his data to the 
Worse-Schimdt solution [Ref. 36], for heat transfer in tubes 
and annular passages. In both the Durao and Ballard studies, 
the experimental data compared relatively well to the predicted 
solutions of MCuen, et al. A curve of the predicted solution 
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of the McCuen, et. al. study is shown in Figure 10. The re¬ 
strictions of that study were as follows: 

1. the velocity profile is fully established. 

2. the fluid transport properties and density are assumed 
constant. 

3. axial heat conduction is negligible. 

4. viscous energy dissipated is negligible. 

5. conditions are invarient with time (steady state). 

All these restrictions have been accounted for in the 

design of the experimental configuration and/or verified in 
the experimental processes. A dimensionless axial length 
coordinate defined as: 



was calculated for 'x' equal to 29.2 centimeters (the length 
of the straight test section). In this equation, ’Pr 1 is 
the Prandtl number of air, taken as 0. 7 , '' is the hydraulic 
diameter (1.239 centimeters), and 'Re ' is the Reynolds number, 
based on the hydraulic diameter. 

Based on the values of x*, it was determined that the 
experiments, for laminar flows, were in the thermal entrance 
region. It was therefore decided to integrate the expression 
of the local Nusselt number as given by McCuen, et. al. [Ref. 
35] : 
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for the small values of x, in order to obtain the following 
expression for the average Nusselt number: 

3 


Nu = - 


4(1 


? Q • -(2/3) 

— ln (1 " 076709T C1 + q h1)x * } 


^WO 


-)x ! 


In this expression, Q . was the heat convected to the fluid 

W 1 

from the unheated wall, and Q was the heat convected to the 

'wo 

fluid from the heated wall. In accordance with the previous 
assumptions, Q . was considered equal to zero, and Q was 
considered equal to Q . , and the average Nusselt number could 

a. 1 l 

then be calculated, and the expression reduces to: 


Nu « In (1 - 


oT6?'0"9 5' 


,-(2/3) 


) 


This value of average Nusselt number was plotted for values 
corresponding to an x* less than 0.042, as this represents the 
maximum value for which the analytic solution for Nusselt 
number is valid. This value of x*, (0.042), for x equal to 
29.2 centimeters, equates to a Reynolds number value of 
approximately 800 , based on hydraulic diameter (or 411 based 
on the channel height). 

The graphical representation shows that the experimental 
results of this study as well as those of Ballard and Durao 
plotted above the analytical solution. This difference can 
probably be explained, at least in part, by the difference in 
geometries used in the studies, and by the limitations inherent 
in any experimental work. The side wall effects and the in¬ 
ability to totally account for all the heat transfer processes 







Comparison of present data with Durao, Ballard and McCuen for Laminar 
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and/or losses could account for the differences of approximately 
nine to ten percent. 

In the comparison of the experimental results obtained in 
this study for those flows not in the laminar regime, the 
comparison and correlation was made to those analytical solu¬ 
tions and experimental results for completely turbulent flows. 
For the most part, the analytical solutions are valid for 
Reynolds number values greater than 10 (equating to a Reynolds 
number of 5128 in this study), and tube flow (or other geo¬ 
metries based on hydraulic diameter). 

A graphical comparison of the results obtained in the 
straight test section runs with the analytical solutions and 
other experimental results is made in Figure 12. The plot 
shows the curve of heat transfer in a straight tube for con¬ 
stant wall temperature, as given by the Dittus-Boelter equations 
[Ref. 38]: 

Nu = 0.023 Re 0,8 Pr 0,4 

and by the empirical relation as given by E.R.G Eckert [Ref. 


39]: 


Nu = 


0.0396 Re 0,75 Pr 


(1 - 


Re 


—7 

T7UZ ) 


This analytical expressions when extrapolated down show some 
agreement with the values obtained at flow rates between forty 
and seventy percent of the rotameter full scale. Unfortunately, 
those values obtained at the higher flow rates showed less of 
a correlation than those mentioned above. These results, 
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obtained for volumetric flow rates higher than seventy percent 
may have been affected by the fact that they were near the top 
end of the rotameter scale, and therefore, may inherently be 
somehwat less accurate than the mid-scale readings. The ther¬ 
mocouple problem discussed previously, may also have contri¬ 
buted to the decrease in correlation. 

When extrapolated down, the experimental results of Kays 
and Leung [Ref. 28], for heat transfer in annular passages, 
with one wall heated at a constant heat rate, and r* equal to 
1.0, (where r* is the annulus radius ratio r^/r ), equivalent 
to parallel plates, appears to correlate well with the present 
results. This extrapolation, was again utilizing a least 
squares straight line correlation, as given in Appendix C. 

Also shown are the analytical solutions of Petuhkov and Popov 
[Ref. 25], for turbulent heat transfer in round tubes given 
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- 

i 
I 

I 

The analytical expression, given by Shibani and Ozisik [Ref. 

26], for turbulent heat transfer between parallel plates: 

I 

Nu. - 12 ♦ 0.03 Re 3 Pr b 

where, J 

a - 0.3 8 - -— 

(3.6 ♦ Pr) 

b * 0.33 ♦ 0.5 e ' 0,6 Pr 

valid for 0.1<Pr <10, and 10^<Re <10 b , is also plotted. 

While no conclusions can be drawn from the data obtained 
because it is not in the range of Reynolds numbers of these 
studies, it is worthy of note, that the present data as extra¬ 
polated, and the analytical and experimental works of others 
can be graphically shown to be close. It should also be noted 
that the investigations of Hsu [Ref. 34], Ashton [Ref. 33, and 
Haynes and Ashton [Ref. 32] all plotted experimental results 
between the Petuhkov and Popov expression and that of Shibani 
and Ozisik, as shown graphically in the Haynes and Ashton 
study above. 

The curved section data is plotted in Figure 13, and 
compared to the experimental work of Brinich and Graham [Ref. 

31] for flow and mass transfer in a curved channel. The 
analytical solution of Petukhov and Popov, and Shibani and 
Ozisik are also plotted as a reference. A least squares 
correlation as given in Appendix C was used to extrapolate 
the data of Brinich and Graham. The accuracy of the data 
points used in this correlation are subject to errors, in that 






the actual values were not given in their study and had to be 
taken from a plot of Stanton number versus Reynolds number as 
given in the reference. 

Irrespective of the relatively small differences in the 
results and the correlations, the trends appear to be conclu¬ 
sive. The heat transfer rates can be seen to increase with 
increasing Reynolds number, and also that the heat transfer 
rate in a curved channel is measurably increased over that 
rate in a curved channel is measurably increased over that 
rate in a straight channel. This increase in heat transfer 
rate, appears to begin with the development of the Taylor- 
Goertler vortices and continues into the transition and tur¬ 
bulent flow regions. This result also appears to have been 
reached by K. V. Dement'eva and I. Z. Aronov [Ref. 41], who 
in 1978 , investigating hydraulic resistance and heat transfer 
in curvilinear channels of rectangular cross-section, 
concluded that there was an enhancement of heat transfer 
intensity in curved rectangular corss-section channels as 
compared to straight pipes and channels. 
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VI. RECOMMENDATIONS 


The amount of experimental work remaining in the study of 
heat transfer in curved channels is significant. Additional 
experiments in the transition region and in the fully turbu¬ 
lent region of fluid flow, should be conducted to refine the 
correlations for average Nusselt number. With this work, it 
would be necessary to design and construct a larger experi¬ 
mental configuration, capable of higher volumetric flow rate 
and greater temperature differences between the heated surface 
and the fluid inlet temperature. Additionally, care should 
be taken to ensure that the temperature of the exit flow, of 
each test section, is accurately measured. This may be 
accomplished by the use of a mixing cup arrangement at the 
exit flow locations. In a channel designed and constructed 
as the one in this study, care must be taken to allow for the 
removal of the mixing cup at the straight section exit location, 
to ensure that turbulence is not induced in the flow by the 
thermocouples prior to the fluid entering the curved section. 

Also in the way of temperature measurement, a thermocouple 
should be installed at the entrance to the rotameter, or pre¬ 
ferably some better device or means of determining the 
volumetric flow rate. This would ensure a more accurate 
measurement of the fluid temperature for an evaluation of the 
fluid properties in the mass flow calculation. 




















= .0200 

The uncertainty in the mass flow rate was then calculated as; 

^ = .0200 

m 

and the Reynolds number uncertainty was 

- 

\ (.0200) 2 + C.01987) 2 + (.00223) 2 + (.01976) 2 = .0345 
d 

Re d = 2042 + 71 

The uncertainty in the Nusselt number was also calculated 

in an identical manner, using the equation: 
tr-T _ Ed 


and the uncertainty as 

) 2 + (|p ) 2 + (S^il ) 2 

Nu » E K air 

It was therefore, necessary to first calculate the uncertainty 
for E. Using the same approach, 

11 a . At 


dh _ J r dQair, 2 . , PL,2 . ,dAT,2 

r ^ 

where AT = T - T„ 
wo B 


■ .0026 and S-I ■ .0048 
The uncertainty in Qair was calculated by: 


Qair * mC pair CT qut - T in ) 
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and 


Qair 
where 


dQair J,dm, 2 . , dC Pair^2 . , d(;T OUT T IN),2 

—*— ~ V + It- > + It -rr- > 

' ' L Pair OUT IN 


m 


and 


— = 0.0200 

m 

d(T 


dC 


Pair 


"Pair 


= .0042, 


OUT ' T IN ^ 


T OUT ' t in 


= .0623 


These calculations result in the following uncertainties 
gfil = .0656 

— = .0658 


h 

dNu 

Nu 


= .0688 


It is therefore, the uncertainty in the difference between 
outlet and inlet temperatures of the flowing air, that is 
ultimately the major source of average Nusselt number uncer¬ 
tainty. This then is the basis for the recommendation to 
better ascertain the exit temperatures for each test section. 
The calculations, as shown above, were for the 651 straight 
section run. Similar calculations for the curved section at 
65% yield the following; 
dRe 


— = .0345 


U 

Re d = 2038 + 70 


dNu 

Nu 


.0563 
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Nu = 9.34 + 0.525. 

The values for the uncertainties for all the variables 
the curved section test run are as follows: 


QUANTITY 

UNCERTAINTY 

a pl 

.0026 

Ac 

.0198 

C P 

.0042 

d 

.0199 

De 

.0398 

FT 

.0526 

Kair 

.0003 

Nu 

.0563 

• 

m 

.0200 

Q 

.0200 

Qair 

.0523 

Re 

. 0345 

t b 

.0073 

t in 

.0029 

t out 

.0106 

T 

W 

.0060 

T -T 

OUT IN 

.0482 

AT 

.0047 

U 

.0022 

P 

.0007 
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APPENDIX B: SAMPLE CALCULATIONS 


A sketch of the control volume for the energy balance, 
for the straight section, indicating the major heat transfer 
components is shown in Figure 14. The sample calculations 
that follow, show the methods and equations that were used 
to obtain these heat transfer components. The heat transfer 
components were required to satisfy the energy balance, and 
to obtain an average Nusselt number as a function Reynolds 
number of the flowing air. 

The curved section energy balance and calculations are 
similar 



FIGURE 14. Energy balance in straight section 
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A. SAMPLE CALCULATIONS DATA 


t room 

= 

21.6°C 


Q 

= 

55% 


V PR 

= 

3.2165 V 


V H 

= 

37.1210 V 


r pr 

= 

2.078 


T 

IN 

= 

0.839 mV = 

21.82° C 

T INS1 

= 

1.324 mV = 

33.34°C 

T INS2 


1.583 mV = 

39.45°C 

Twol 

= 

1.840 mV 


Two 2 

= 

1.840 mV 


Two3 

= 

1.987 mV 


Two4 

= 

2.001 mV 

Two = 

Two 5 

= 

2.112 mV 


Two6 

= 

2.125 mV 


Two 7 

3 

2.245 mV 


Two8 

= 

2.241 mV 


Twil 

= 

0.964 mV 


Twi2 


0.967 mV 


Twi3 

= 

1.012 mV 


Twi4 

3 

1.008 mV 

Twi = 

Twi5 

3 

1.056 mV 


Twi6 

3 

1.047 mV 


Twi7 

3 

1.101 mV 


Twi8 

* 

1.100 mV 


t out 

3 

1.124 mV 

= 27 . 97°C 

t exit 

— 

1.076 mV 

= 26 . 91°C 

a pl 

= 

0.0742 m 2 



50.09°C 


25.76°C 


= 80.44° 
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A 

c 

= 

0.0016 m 

d 

= 

0.00635 m 

K air 

= 

0.02662 w/m°C 

r 

^Pair 

u air 

= 

1.0057 KJ/Kg°C 
1.983 x 10' 5 Kg/m 

k ins 

= 

4.18 x 10 W/m°C 

ax ins 

s 

0.00635 m 

£ 

wo 

= 

0.70 

c 

"Wl 

3 

0.40 

F 

wo-wi 

a 

: 

1.0 

5.669 x 10" 8 W/m 2 


TEMPERATURE CALCULATIONS 

1. Temperature Difference in Insulation (AT^g) 

AT INS = T INS2 ’ T INS1 * 39,45 - 33.34 = 6.11°C 

2. Bulk Temperature (T g ) 

Tb , T IN * .. T OUI , 21 -82 * 27.97 . 24 _ 90 . c 

3. Mean Temperature Difference (AT) 

AT = T - T- - 50.09 - 24.90 = 25.19°C 
wo B 


POWER CALCULATIONS 


1 . 


2 . 


Power Supplied (Q^) 

« . V PR V H _ (3.2165)(37.121) 

Q P — C2TUTF)- 


57.46 W 


Heat Lost Through Outer Plate (Q^ ) 

q - AT INS = _ 6.11 __ 

Lo AX INS /K INS' A PL (0.00635)/((4.1 8x10 _Z )(0.0742)) 


=■ 2.98 W 
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ft 


3. 


Heat Radiated (Q ) 
a. Radiation Resistance (R R ) 


1 -e 


n = WO 

R A PL £ wo 


, 1 -e . 

+ 1 + wi 

+ A F T + A £ 
A PLwo-wi A PL wi 


+ _i_ - 1] = 2^929 
A PL e wo £ wi A PL 


b. Heat Radiated (Qr) 
a (T 4 -T. . 4 ) 

WO Wl J 


Qr = 


'K 


R, 


C5.669x10" 8 )(0.0742)(323.25 4 - 
2.929 


= 4.214 W 

Heat Convected to Air (Q . ) 

_ all 

a. Corrected Volumetric Flow Rate (q) 

t fxtt (°r) 

q = q/0 = q * rnrr 

Standard 1 


(0.793) (0.55) 
60 


540.13 

529.69 


b. 


= 0.0073m J /sec 
Density (p ) 


298.92 4 ) 


P 


101325 






c. Mass Flow Rate (m) 


m = qp = (0.0073)(1.177) = 0.0086Kg/sec 
d. Heat Convected to Air (Q . ) 

all 

Q air = m C Pair (T OUT ' T IN ) 

= (0.0086)(1.0057)(27.97-21.82) 

= 0.05319 KJ/sec 

= 53.19 W 


D. AVERAGE HEAT TRANSFER COEFFICIENT (K) 

tr _ Q air _ 53.19 

h ■ a pl AT “ Co'on2) (25.13") 

= 28.46 W/M 2o C 


E. AVERAGE NUSSELT NUMBER (Nu) 


Nu 


Ed _ (28.46)(0.00635) _ 
Kair (0. 02662) 


6.789 


F. REYNOLDS NUMBER (Re) 

Re = m d = (0.0086)(0,00635) 
d y air A c (1.983xl0' 5 )(0.0016) 


1721.19 




APPENDIX C: CORRELATIONS 


All correlations were obtained using the method of least 

squares (first degree polynomial), as outlined by C. F. Gerald 

[Ref. 43]. In this method the values of Reynolds number (x) 

and Nusselt number (Y) were first converted to the appropriate 

natural logrithemic value, after which the summations of x., 

2 

x i ’ Y i> anc * x i Y i were calculated. The values of these quanti¬ 
ties were then placed in a matrix to solve the simultaneous 
equations: 

a o N * ajSXj - SYj 

a o Ix i * ’ Ix i Y i 

where N is the number of data points, and the matrix is 


N 

s *i’ 


a 

0 


’* Y 1 ‘ 

2X i 



a l 

S 

5X 1 Y 1 


After attaining the values of a Q and a^, the equation of 

the line became: 

Ln Nu * a + a.LnRe 
o 1 

or 

_ a a, 

Nu * e Re 

The quantities for the various studies are given in Table 






The standard deviation for the present study curved section 


correlation was calculated as: 



where: is the actual value of average Nusselt number ob¬ 

tained in the test runs 


y^ is the value of average Nusselt number calculated 
from the correlation equation 

N is the number of data points used in the correlation 
n is the degree polynomial used. 

The result of the above method was 
a 2 = .02426 
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TABLE III: CORRELATIONS 



V 



! 










H 
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H 

O 
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HH 
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HH 
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HH 
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2 
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2 

^3 

so 


HH 
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03 
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(2Q 
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2 
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2 
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OS 

tu 

►j 
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LO 
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vO 
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LO 
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'O 
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